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ABSTRACT 

Aims. We present our analysis of a local ( z = 0.04 — 0.2) sample of 31 galaxy clusters with the aim of measuring 
the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set 
constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. 

Methods. We have exploited the large field-of-view and low instrumental background of ROSAT /PSPC to trace the 
density of the intracluster gas out to the virial radius. We stacked the density profiles to detect a signal beyond r2oo and 
measured the typical density and scatter in cluster outskirts. We also computed the azimuthal scatter of the profiles 
with respect to the mean value to look for deviations from spherical symmetry. Finally, we compared our average density 
and scatter profiles with the results of numerical simulations. 

Results. As opposed to some recent Suzaku results* and confirming previous evidence from ROSAT and Chandra , we 
observe a steepening of the density profiles beyond ~ r^oo. Comparing our density profiles with simulations, we find that 
non-radiative runs predict density profiles that are too steep, whereas runs including additional physics and/or treating 
gas clumping agree better with the observed gas distribution. We report high-confidence detection of a systematic 
difference between cool-core and non cool-core clusters beyond ~ 0.3r2oo, which we explain by a different distribution 
of the gas in the two classes. Beyond ~ 7 * 500 , galaxy clusters deviate significantly from spherical symmetry, with only 
small differences between relaxed and disturbed systems. We find good agreement between the observed and predicted 
scatter profiles, but only when the 1% densest clumps are filtered out in the ENZ0 simulations. 

Conclusions. Comparing our results with numerical simulations, wc find that non-radiative simulations fail to reproduce 
the gas distribution, even well outside cluster cores. Although their general behavior agrees more closely with the 
observations, simulations including cooling and star formation convert a large amount of gas into stars, which results in 
a low gas fraction with respect to the observations. Consequently, a detailed treatment of gas cooling, star formation, 
AGN feedback, and consideration of gas clumping is required to construct realistic models of the outer regions of 
clusters. 

Key words. X-rays: galaxies: clusters - Galaxies: clusters: general - Galaxies: clusters: intracluster medium 


1. Introduction 

The outskirts of galaxy clusters are the regions where 
the transition between the virialized gas of clusters 
and the accreting matter from large-scale structure oc- 
curs and where the current activity of structure for- 
mation takes place. Around the virial radius, the as- 
sumption of hydrostatic equilibrium, which is a neces- 
sary assumption for reconstructing cluster masses from 
X-ray measurements, might not be valid any more (e.g., 
Evrard et al. 1996), which could introduce biases to X- 
ray mass proxies (Rasia et al. 2004; Piffaretti & Valdarnini 
2008; Nagai et al. 2007b; Lau et al. 2009; Meneghetti et al. 
2010; Fabjan et al. 2011). As a result, the characterization 


of the X-ray emitting gas in the outer regions of galaxy clus- 
ters is important for mapping the gas throughout the entire 
cluster volume, studying the formation processes currently 
at work in the Universe, and performing accurate mass es- 
timates for cosmological purposes (e.g., Allen et al. 2011). 


Because of the low surface brightness of the X-ray 
emitting gas and the extended nature of the sources, 
measuring the state of the intracluster gas around the 
virial radius is challenging (Ettori &; Molendi 2011). 
Recently, the Suzaku satellite has achieved a break- 
through in this domain, performing measurements of 
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cluster temperatures out to 7 * 200 1 (Reiprich et al. 2009; 
Bautz et al. 2009; Kawaharada et al. 2010; Hoshino et al. 
2010; Simionescu et al. 2011; Akamatsu et al. 2011; 
Humphrey et al. 2012), even in one case beyond 7*200 
(George et al. 2009), although this detection is likely 
hampered by systematic effects (Eckert et al. 2011a). 
Interestingly, some of the Suzaku results indicate very 
steep temperature profiles and shallow density profiles 
in cluster outskirts, at variance with the results from 
XMM-Newton (Pratt et al. 2007; Leccardi &; Molendi 
2008, Snowden et al. 2008; Croston et al. 2008), Chandra 
(Vikhlinin et al. 2006; Ettori &; Balestra 2009), ROSAT 
(Vikhlinin et al. 1999; Neumann 2005), and with the 
results from numerical simulations (Roncarelli et al. 2006; 
Tozzi & Norman 2001; Nagai & Lau 2011). Thus, the 
behavior of the gas in cluster outskirts is still the subject 
of debate. Throughout paper, we refer to cluster outskirts 
as the region with r > 7 * 500 - 

Thanks to its large field of view (FOV, ~ 2 deg 2 ) 
and low instrumental background, ROSAT/ PSPC is to the 
present day the most sensitive instrument for low surface- 
brightness emission. Its ability to detect cluster emission at 
large radii has been demonstrated by Vikhlinin et al. (1999) 
and Neumann (2005) (hereafter, V99 and N05). Because 
of the large FOV, it can perform simultaneous local back- 
ground measurements, so it is less affected than Suzaku by 
systematic uncertainties. Its main limitation, however, is 
the r estricted band pass and poor spectral resolution, which 
makes it impossible to measure cluster temperatures. 

This paper presents the analysis of a sample of 31 
galaxy clusters observed with ROSAT/ PSPC, with the aim 
of characterizing the cluster emission at large radii and 
comparing the results with three different Sets of numer- 
ical simulations (Roncarelli et al, 2006; Nagai & Lau 2011; 
Vazza et al. 2010). The paper is organized as follows. In 
Sect. 2, we describe our cluster sample and the available 
data. V/e present our data analysis technique in Sect. 3 and 
report our results in Sect, 4. We compare our results with 
numerical simulations in Sect. 5 and discuss them in Sect. 
6 . 

Throughout the paper, we assume a ACDM cosmology 
with n m = 0.3, n A = 0.7, D b = 0.047, and H 0 = 70 km s" 1 
Mpc -1 . 

2. The sample 

We selected objects in the redshift range 0.04 — 0.2, such 
that r 2 oo is easily contained within the FOV of the instru- 
ment and is large enough to allow for an adequate sam- 
pling of the density profile. We restricted ourselves to ob- 
servations with enough statistics to constrain the emission 
around the virial radius. Our final sample comprises 31 clus- 
ters in the temperature range 2.5-9 keV, with the addition 
of A2163 ( kT ~ 18 keV). Among our sample, we classified 
14 clusters as cool core (CC) following the classification 
of Cavagnolo et al. (2009) (i.e. they exhibit a central en- 
tropy Ko < 30 keV cm 2 ), and 17 as non cool core (NCC, 
Kq > 30 keV cm 2 ). We recall that CC clusters exhibit a 
relaxed morphology, a high central density and a temper- 
ature decrement in the central regions, while NCCs trace 
dynamically-disturbed clusters with irregular morphologies 

1 We define ta as the radius within which M(< ta)/^t\ - 

A p cr ii 


and flat temperature and density profiles in their cores (e.g., 
Sanderson et al. 2009; Hudson et al. 2010). 

Our sample of clusters, together with the log of the 
available data and some important quantities, is shown in 
Table 3. In Fig. 1 we plot the distribution of temperature 
(left hand panel) and central entropy (right hand panel) 
for our sample. It should be noted that the sample was 
selected based on the quality of the existing observations 
and might be subject to selection effects. However, for the 
purpose of this work we did not require that the sample 
be representative or complete, since we are interested in 
characterizing cluster outskirts, which exhibit a high level 
of self-similarity. 


3. Data analysis 

3.1. Data reduction 

We used the ROSAT Extended Source Analysis Software 
(Snowden et al. 1994) for data reduction. We filtered out 
time periods when the master veto count rate exceeds 220 
cts/sec (using valid-times), and extracted light curves for 
the whole observation using rate_pspc. We used the ao 
executable to model the atmospheric column density for 
the scattering of solar X-rays, and fit the light curves in 
each energy band to get the relative contributions of the 
scattered solar X-rays (SSX) and of the long-term enhance- 
ments (LTE), using the rate Jit executable. 

We then extracted event images in each energy band 
and the corresponding effective exposure maps, taking vi- 
gnetting effects into account. We computed the contri- 
bution of the various background components, the LTE 
(lte_pspc), the particle background (cast_part), and the 
SSX (cast_ssx), and combined them to get a map of all 
the non-cosmic background components. 

3.2. Surface-brightness profiles 

The point-spread function (PSF) of ROSAT /PSPC 
strongly depends on angle, and ranges from ~ 15 arcsec on- 
axis to 2 aremin in the outer parts of the FOV. Thus, the 
sensitivity of the instrument to point sources is higher on- 
axis, and a larger fraction of the cosmic X-ray background 
(CXB) is resolved. Consequently, when detecting sources in 
the image it is important to use a constant flux threshold, 
such that the same fraction of the CXB is resolved over 
the entire FOV and the value measured in the source-free 
regions can be used to subtract the background. We detect 
point sources using the program detect with a minimum 
count rate of 0.003 cts/sec in the R3-7 band (~ 3 x 10~ 14 
ergs cm' 2 s -1 in the 0. 5-2.0 keV band) to resolve the same 
fraction of the CXB over the FOV, and mask the corre- 
sponding areas. To compute surface-brightness profiles, we 
extract count profiles from the event images in the R3-7 
band (0.42-2.01 keV) with 30 arcsec bins centered on the 
surface-brightness peak, out to the radius of 50 aremin. 
We divide each pixel by its corresponding exposure to ac- 
count for the vignetting effects, following the procedure of 
Eckert et al. (2011b) 2 . We perform the same operation for 
the background map and subtract the non-cosmic back- 
ground profile in each bin. 


2 http:/ /www.iasf- milano.inaf.it / ~eckert/newsite/Proffit.html 
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We tested this procedure on four different blank fields 
to estimate the accuracy in our determination of the CXB. 
We extracted the surface-brightness profile for the four ob- 
servations from the center of the FOV, grouped the bins 
to ensure a minimum of 100 counts per bin, and fitted 
the resulting profiles with a constant (see Fig. 2). While 
the agreement is qualitatively good, significant deviations 
from the model are found, leading to an excess scatter of 
~ 6%, which we used as an estimate of the systematic un- 
certainties in measuring the CXB. This value encompasses 
both the cosmic variance and the true systematic uncertain- 
ties, e.g., in the vignetting correction or determination of 
the particle background. The higher level of scatter in the 
central regions is explained by the small area of the cor- 
responding annuli, which implies a large cosmic variance 
likely due to discrete sources with fluxes just under our ex- 
clusion threshold. Since in most cases the value of 7*200 is 
larger than 15 arcmin, our systematic error of 6% is a con- 
servative estimate of the level of systematic uncertainties 
at the virial radius. 

For each cluster, we then use temperature pro- 
files from the literature ( XMM-Newton , Snowden et al. 
(2008); Chandra , Cavagnolo et al. (2009); BeppoSAX, 
De Grandi & Molendi (2002)) to estimate the virial tem- 
perature of the cluster. We approximated T V { r as the mean 
temperature in the 200-500 kpc region, i.e. excluding the 
cool core and the temperature decline in the outskirts 
(Leccardi &; Molendi 2008). Using this estimate of T viri we 
computed the value of r2oo from the scaling relations of 
Arnaud et al. (2005). One might argue that the scaling re- 
lations of Arnaud et al. (2005) were computed using the 
mean temperature in the 0.1 -6.67*200 region, which in most 
cases extends beyond the available temperature profiles. 
Using the mean temperature profiles of Leccardi & Molendi 
(2008), we computed the mean temperature extracted in 
the 0.1 - 0.5r2oo and 200-500 kpc regions. In the temper- 
ature range of our sample, we found that the results differ 
at most by 2%, so our values of r 2 oo are unbiased. We then 
used the source- free region of the observation (r > 1.37*2oo) 
to fit the surface-brightness profile with a constant and get 
the cosmic background level for the observation, with the 
exception of the Triangulum Australis cluster, for which we 
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Fig. 2. Surface-brightness profiles for 4 blank-field PSPC 
observations from the center of the FOV, fitted with a con- 
stant. The dashed line shows the vignetting correction curve 
for comparison, in arbitrary units; the bump at ~ 22 ar- 
cmin is caused by the support structure. The bottom panel 
shows the ratio between data and model. 


used the range r > l.lr2oo because of the high value of 7*200 

37 arcmin). 

After having estimated the sky background for our ob- 
servation, we again extracted the surface-brightness profile 
in the radial range 0 - 1.3r 2 oo with logarithmic bin size. 
The best-fit value for the CXB was subtracted from the 
profile and its error was added in quadrature to each bin. 
The systematic error of 6% on the CXB was also added in 
quadrature to account for the cosmic variance and system- 
atic uncertainties. For comparison, we note that in most 
cases the statistical uncertainties in the profiles are on the 
order of 10% of the CXB value around r 2 00 • 
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3.3. Density profiles 

To compute the density profiles, we first rebinned our 
background-subtracted surface- brightness profiles to ensure 
a minimum of 200 counts per bin and a detection signifi- 
cance of at least 3<r , to reach sufficient statistics in each bin. 
We used the procedure of Kriss et al. (1983) to deproject 
the observed profiles, and the PSPC response to convert 
the observed count rates into emission measure, through 
the normalization of the MEKAL model (see Eckert et al. 
2011a, for details), 


Norm — 


1Q~ 14 

47r[d A (l + z)] 2 


I 


n e nf{dV , 


(i) 


which is proportional to the emission measure. We assumed 
that the spectrum of our sources is described by an ab- 
sorbed MEKAL model with Nh fixed to the 21cm value 
(Kalberla et al. 2005) and abundance fixed to O.SZq. We 
used temperature profiles from the literature (see Table 3) 
and interpolated them onto the same grid as the SB pro- 
files. The resulting model was then folded with the PSPC 
response, and the conversion from PSPC R3-7 count rate 
to emission measure was inferred. Beyond the limit of the 
temperature profiles, the temperature of the outermost an- 
nulus was used. We note that the conversion from PSPC 
count rate to emission measure is highly insensitive to the 
temperature; between 2 and 8 keV, the conversion factor 
changes at most by 4%. Once converted into the MEKAL 
normalization, we inferred the density profiles, assuming 
spherical symmetry and constant density into each shell. 

The error bars on the density profiles were estimated us-, 
ing a Monte Carlo approach. In each case, we generated 10 4 
realizations of the surface-brightness profile using Poisson 
statistics, and performed the geometrical deprojection fol- 
lowing the method described above. The ler error bars were 
then estimated bv computing the root-mean square devia- 
tion (RMS) of our 10 4 realizations of the density profile in 
each density bin. 


3.4. Azimuthal scatter profiles 

For the purpose of this work, we are also interested in the 
deviations in the X-ray emission from spherical symmetry. 
We divide our images into N azimuthal sectors with con- 
stant opening angle, and compute the surface-brightness 
profiles in each sector individually. We then compute the 
scatter of the various sectors with respect to the mean 
profile, following the definition introduced by Vazza et al. 
(2011b), 


2 1 ^ (SBj - (SB)l 

N £ (SB)* ’ 


( 2 ) 


where (SB) is the mean surface-brightness and SBi,i = 
1..N denotes the surface-brightness computed in the vari- 
ous sectors. It must be noted that the statistical fluctua- 
tions of the SB between the different sectors introduce a 
certain level of scatter in Eq. 2, which must be taken into 
account for determining the level of intrinsic scatter. We 
used two different methods to disentangle between statis- 
tical and intrinsic scatter. In the first case, we computed 
the level of statistical scatter independently and subtracted 


it from Eq. 2. In the second case, we used a maximum- 
likelihood estimator to determine the intrinsic scatter and 
its uncertainties. The two methods gave consistent results 
and are described in detail in Appendix A. For the remain- 
der of the paper, we refer to the results obtained using the 
direct method (see Sect. A.l). 

In our analysis, we group the bins of the total surface- 
brightness profiles to reach a minimum of 8<r per bin to en- 
sure adequate statistics in the scatter measurements, and 
then divide our images into 1 2 sectors with an opening of 
30°. The result of this analysis is a radial profile describing 
the intrinsic azimuthal scatter of the X-ray surface bright- 
ness, in percent. 

It must be noted that the method presented here is sen- 
sitive to all kinds of deviations from spherical symmetry, 
whether it is induced by the asymmetry of the large-scale 
structure (e.g,, filaments), by gas clumping or by ellipticity 
The cause of the observed asymmetry cannot be determined 
from the azimuthal scatter alone. 


4. Results 

4.1. Emission measure and density profiles 

In Fig. 3 we show the scaled emission measure profiles (left, 
following Eq. 1) and the deprojected densit\ profiles (right) 
for the 31 clusters in our sample. A self-similar scaling 
was applied to the emission-measure profiles (Arnaud et al. 
2002); i.e., each profile was rescaled by the quantity 


( 3 ) 

The density profiles were rescaled by E 2 (z) = H m (l -f 
z) 3 4- 11a following their expected evolution with redshift 
(Croston et al. 2008). As already noted by several authors 
(e.g., Vikhlinin et al. 1999; Neumann 2005; Croston et al. 
2008; Leccardi et al. 2010), the profiles show a remarkable 
level of self-similarity outside of the core (r > Q.2r2oo)- 
On the other hand, the large scatter observed in the cen- 
tral regions reflects the distinction of the cluster popu- 
lation into CCs, showing a prominent surface-brightness 
peak, and NCCs, which exhibit a flat surface brightness pro- 
file in their cores, as expected from the standard model 
(Cavaliere &; Fusco-Femiano 1976), 


SB(r) - SB 0 



( 4 ) 


In the radial range 0.2 — 0.7r2oo, the scatter of the den- 
sity profiles is 10%-20%, in excellent agreement with the 
Chandra (Vikhlinin et al. 2006) and X MM- Newton re- 
sults (Croston et al. 2008). However, Croston et al. (2008) 
needed to rescale the profiles by T” 1 / 2 to account for the 
lower gas fraction in low-mass objects. In our case, perform- 
ing such a scaling does not reduce the scatter of the profiles 
further. This is probably explained by the relatively narrow 
temperature range spanned in our sample (all but one ob- 
jects have a temperature higher than 3 keV), such that the 
clusters in our sample should show little dependence on gas 
fraction. 
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Fig. 3. Scaled emission measure (left, in units of cm 6 Mpc) 
(see Table 3). 

4.2. Stacked emission-measure profiles 

To compute the mean profile of our sample, we interpolated 
each profile following a predefined binning in units of 7*200 
common to all clusters and performed a weighted mean to 
compute stacked profiles. The errors on the interpolated 
points were propagated to the stacked profiles. We also di- 
vided our sample into the two classes (CC and NCC) to 
look for differences between them. 

In Fig. 4 we show the stacked emission-measure (EM) 
profile for the entire sample compared to the profiles 
stacked for the two populations separately (see also 
Appendix C). Interestingly, we note a clear distinction be- 
tween the two classes in cluster outskirts (see the bottom 
panel of the figure). Namely, beyond ~ 0. 37 * 200 > NCC pro- 
files systematically exceed CCs. A similar effect has re- 
cently been noted by Maughan et al. (2011), who found 
a crossing of the average density profiles at a similar ra- 
dius, and also at a lower statistical significance in the works 
of Amaud et al. (2010) and Pratt et al. (2010). We stress 
that this effect is really a difference between the two classes; 
i.e. it is not introduced by a biased distribution of another 
quantity (such as temperature or redshift). Indeed, group- 
ing the profiles according to the temperature or the redshift 
did not show any particular behavior, which indicates that 
we are really finding an intrinsic difference between the CC 
and NCC classes. This result could follow from a different 
distribution of the gas in the two populations or from a 
higher clumping factor in disturbed objects (see Sect. 6). 

Alternatively, the observed difference could be explained 
by an inaccurate determination of 7*200 for NCC clusters. 
Indeed, the scaling relations of Arnaud et al. (2005) were 
computed under the assumption of hydrostatic equilibrium, 
which is fulfilled better in CC clusters. This explanation 
is, however, unlikely. Indeed, to recover self-similarity, our 
value of 7*200 should have been systematically underesti- 
mated by ~ 10% for NCCs, i.e. since 7*200 oc the virial 
temperature of the NCC clusters should have been under- 
estimated by more than 20%. From mock Chandra observa- 
tions of a sample of simulated galaxy clusters, Nagai et al. 
(2007b) have determined that the spectroscopic tempera- 
tures of unrelaxed clusters differs from that of relaxed clus- 
ters by ~ 5%, which is not enough to explain the observed 



and density profiles (right) for the 31 dusters of our sample 



Fig. 4. Stacked emission measure profile (in units of cm 6 
Mpc) for the entire sample (black), and the two populations 
individually (CC, red; NCC, blue). See also Appendix C. 
The bottom panel shows the ratio between the CC and 
NCC populations. 


difference. It is therefore unlikely that such a large error on 
the virial temperature would be made. 

We fit the mean scaled emission-measure profiles from 
Fig. 4 with the standard £-modei (Eq. 4), adding a second 
P component in the case of the CC clusters to take the 
cool core into account. The (double) ft model gives a good 
representation of the data in the radial range 0 - 0.7t*2oo 
(~ 7 * 500)7 but significantly exceeds the observed profiles 
abo-e this radius, in agreement with the results of V99, 
N05, and Ettori h Balestra (2009). For CC clusters, the 
best-fit model gives p = 0.717± 0.005, while for NCC clus- 
ters we find p = 0.677 ± 0.002. Fitting the radial profiles in 
the range 0.65-1.3t*2ooi we observe a significant steepening, 
with a slope P = 0.963±0.054for CCs and P = 0.822±0.029 
for NCCs. As explained above, the slope of the NCC profile 
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Fig. 5. Average proton density profile for the entire sample. 
The dashed lines indicate the positive and negative scatter 
of the profiles around the mean value. 


is flatter than that of the CC profile beyond 7*500 * The fits 
of the profiles in various radial ranges are reported in Table 
2 to quantify the steepening. 

Given the limited number of objects in our sample, we 
have to verify that this result is not a chance realization. 
We fit all the emission-measure profiles at r > 0.3t*2oo with 
a 0 profile, fixing the value of 0 to 0.7 and r c to 0.12r2oo> 
and extracted the best-fit normalization for all profiles. We 
then sorted the normalization values into the CC and NCC 
classes, and performed a Kolmogorov- Smirnov test to de- 
termine the probability that they originate in the same par- 
ent distribution. Using this procedure, we found that the 
chance probability for this result is very low, 6xl0 -7 . 
Therefore, we can conclude with good confidence that we 
are indeed finding an intrinsic difference between the two 
classes. 

4.3. Stacked density profiles 

We stacked the density profiles shown in the right hand 
panel of Fig. 3 following the same method as for the EM 
profiles. From the different profiles, we also computed the 
scatter of the profiles around the mean value, following a 
method similar to the one presented in Sect. 3.4 for the 
azimuthal scatter. The statistical scatter was subtracted 
from the total scatter using the same technique. In Fig. 
5 we show the average density profile of our clusters to- 
gether with the scatter of the individual profiles around 
the mean value (see also Table 1). At 7 * 200 , the mean den- 
sity is ri 2 oo = (3-8 ± 0,4) x 10 ~ 5 E 2 (z) cm” 3 , with 25% 
scatter. For comparison, it is interesting to note that the 
density of PKS 0745-191 claimed in the Suzaku analysis of 
George et al. (2009) at 7*200 deviates from our mean value 
by more than 5a, which casts even more doubt on this mea- 
surement (Eckert et al. 2011a). 

As for the EM, we also extracted mean density profiles 
individually for the two classes of clusters in our sample. 
The same behavior is observed at large radii; i.e., the den- 
sity of NCC clusters is systematically higher (by ~15%) 
than that of CCs above r ~ 0.37*2oo- A global steepening of 
the density profiles is also observed beyond ~ r 50 o- 


Table 1. Mean emission-measure and density profiles com- 
puted from our sample. 


Rin 

Rout 


ScEM 



n H B(z) * 

a 

0 

0.03 

(1.78 

± 

0.01) 

10 

-5 

11.447 

' d: 

: 0.033 

58 

0.03 

0.06 

(1.23 

± 

0.00) 

10 

-5 

6.325 

zb 

0.018 

46 

0.06 

0.09 

(7.34 

± 

0.03) 

10‘ 

-6 

3.446 

± 

0.012 

36 

0.09 

0.12 

(5.13 

± 

0.02) 

10 

-6 

2.222 

± 

0.010 

26 

0.12 

0.15 

(3.49 

± 

0.01) 

10 

-6 

1.599 

± 

0.009 

21 

0.15 

0.18 

(2.44 

± 

0.01) 

10 

-6 

1.191 

± 

0.008 

17 

0.18 

0.21 

(1.65 

+ 

0.01) 

10 

-6 

0.923 

± 

0.007 

17 

0.21 

0.24 

(1.24 

± 

0.01) 

10* 

-6 

0,731 

it 

0.007 

13 

0.24 

0.27 

(9.66 

± 

0.06) 

10' 

-7 

0.606 

± 

0.006 

12 

0.27 

0.30 

(7.19 

± 

0.05) 

10" 

-7 

0.506 

± 

0.006 

15 

0.30 

0.33 

(5.50 

± 

0.04) 

10" 

-7 

0.422 

± 

0.005 

12 

0.33 

0.37 

(4.20 

± 

0.04) 

10 

-7 

0.360 

± 

0.005 

15 

0.37 

0.42 

(3.08 

± 

0.03) 

10“ 

-7 

0.289 

± 

0.005 

12 

0.42 

0.47 

(2.11 

± 

0.02) 

10- 

-7 

0.227 

± 

0.004 

12 

0.47 

0.52 

(1.53 

± 

0.02) 

10- 

-7 

0.193 

± 

0.004 

18 

0.52 

0.59 

(1.05 

± 

0.02) 

10 

-7 

0.143 

± 

0.004 

16 

0.59 

0.66 

(7.16 


0.15) 

10 

-8 

0.121 

± 

0.004 

25 

0.66 

0.74 

(5.12 

± 

0.14) 

10' 

-8 

0.092 

± 

0.003 

10 

0.74 

0.83 

(3.36 

± 

0.12) 

10" 

-8 

0.072 

± 

0.003 

34 

0.83 

0.93 

(1.97 

± 

0.12) 

10 

-8 

0.059 

± 

0.002 

17 

0.93 

1.05 

(1.06 

i 

0.11) 

10 

-8 

0.039 

± 

0.002 

11 

1.05 

1.17 

(6.33 


1.01) 

10 

-9 

0.028 

± 

0.002 

22 


Note: Column description. 1 and 2: Inner and outer bin radii in 
units of 7*2001 3: Emission measure rescaled by A ssc in units of 
cm -6 Mpc; 4: Average proton density in units of 10“ 3 cm~ 3 ; 
5: Scatter of the various profiles relative to the mean value in 
percent. 

Our density profiles are in good agreement with the re- 
sults of V99. However, while V99 estimated the density 
from j#-model fitting, we performed a geometrical depro- 
jection of the data using temperature profiles to infer the 
mean density profile. This method has the advantage of not 
depending on any model. 


4.4. Gas mass 

We computed the gas mass from our deprojected density 
profiles and stacked them in the same way as described 
above. In the self-similar model, the gas mass is expected 
to follow the relation M oc T 3 / 2 (e.g., Bryan & Norman 
1998). However, observational works indicate that the ac- 
tual Mg a s — T relation is steeper than the expected self- 
similar scaling (Neumann k Arnaud 2001; Arnaud et al. 
2007; Croston et al. 2008) because of the lower gas fraction 
in groups and poor clusters. For this work, we use the rela- 
tion determined from the REXCESS sample (Croston et al. 
2008) to rescale our gas mass profiles, 

/ U'T' \ 1986 

As above, we divided the sample into CC and NCC 
classes, and stacked the two classes individually. In Fig. 6 
we show the mean gas mass profiles for CC and NCC clus- 
ters. As expected, CCs have a higher gas mass in their inner 
regions, since their central densities are higher. More inter- 
estingly, we see that the two profiles converge in cluster 
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**200 

Fig. 6. Enclosed gas mass profiles for CC (red) and NCC 
systems (blue). The data were rescaled by E(z)kT~ 1986 as 
observed in the REXCESS sample (Croston et al. 2008). 

outskirts, and exhibit a gas mass around the virial radius 
that is consistent within the error bars. At 7 * 200 , the univer- 
sal gas mass is 

/ wp \ 

M 90S , 200 = (2.41±0.05) x l0 li E(z)~ 1 f J M 0 , (6) 

with a scatter of 17% around the mean value. This result 
follows from the higher density measured in average beyond 
~ 0 . 07*200 in NCC clusters and the steeper slope of CC pro- 
files in the outskirts (see Sect. 4.2). The lower density of 
CC clusters in the outer regions compensates for the well- 
known excess observed in the cores, such that the total gas 
mass contained within the dark-matter halo follows a uni- 
versal relation. We also estimated the average gas fraction 
by computing the expected value of M 200 using the scaling 
relations of Arnaud et al. (2005). For our sample, we find 
a mean gas fraction within 7*200 of 

/ kT \ 0,486 

/ gas, 200 = (0.15 ± 0.01) YoteV) ' ( 7 ) 

in good agreement with previous works (e.g., 
Vikhlinin et al, 2006; McCarthy et al. 2007), which 
for the most massive objects corresponds to ^ 89% of the 
cosmic baryon fraction (Jarosik et al. 2011). 

4.5. Azimuthal scatter 

Following the method described in Sect. 3.4, we computed 
the azimuthal scatter of the surface-brightness profiles for 
all the clusters in our sample, and rescaled the scatter pro- 
files by our estimated value of 7*200- We then stacked the 
profiles using the same procedure as described above and 
computed the mean azimuthal scatter. We recall that since 
the surface brightness depends on n the variations in den- 
sity are less important than the ones computed here. 

In Fig. 7 we plot the average scatter profile, compared 
to the mean value for CC and NCC clusters. The increase 
in the innermost bin is an artifact introduced by the small 
number of pixels in the center of the images, so it should 
be neglected. At small radii (r < 0.5r2oo) we find a clear 



Fig. 7. Stacked azimuthal scatter in surface-brightness for 
the entire cluster sample (black). The red and blue data 
represent the mean profile extracted from CC and NCC 
clusters, respectively. 


difference between CC and NCC clusters, which is eas- 
ily explained by the more disturbed morphology of the 
latter. In this radial range, CC profiles exhibit a scat- 
ter of 20-30%, which corresponds to density variations on 
the order of 10%, in good agreement with the value pre- 
dicted by Vazza et al. (2011b) from numerical simulations. 
Conversely, beyond r ~ 7 * 500 ; the profiles for CC 'and NCC 
clusters are similar, and indicate a high scatter value (60- 
80%). 

We investigated whether any systematic effect could af- 
fect our result in cluster outskirts, where the background 
is dominating with respect to the source. Indeed, in such 
conditions, the total scatter is dominated by the statistical 
scatter. In case the mean level of systematic uncertainties in 
the CXB reconstruction exceeds our adopted value of 6%, 
Eq. A. 2 immediately implies that the intrinsic scatter would 
be overestimated. The presence of both intrinsic and star 
tistical scatter could also introduce some covariance term, 
which is not taken into account in Eq. A. 2. To test this hy- 
pothesis, we ran a set of simulations including source and 
background, where we introduced a given level of intrinsic 
scatter for the source and a systematic error in addition 
to the Poisson statistics for the background. We then com- 
puted the intrinsic level of scatter following Eq. A. 2. Our 
simulations indicate that, even when increasing the level of 
systematic uncertainties to 12% of the CXB value, a signif- 
icant bias in the measured scatter only appears when the 
source-to-background ratio is close to the systematic uncer- 
tainties. Since, by construction, we never detect any signal 
when the source is less than ~15% of the CXB value, our 
results are unaffected by these effects, and we can conclude 
with good confidence that the high level of scatter mea- 
sured beyond ~ 7*500 is an intrinsic property of our cluster 
sample. 

In addition, we also tested whether the scatter for the 
two populations in the outermost regions could be affected 
by small-number statistics or driven by some particular ob- 
jects. Of the 31 objects in our sample, a measurement of the 
scatter at 7*200 could be obtained for 23 of them (12 NCC 
and 11 CC). We used a jackknife method to test whether a 
single object dominates the results for any of the two pop- 
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illations; i.e., we randomly exclude one or two profiles from 
the sample, recompute the mean profiles, and examine the 
distribution of the mean values. In both cases, the distri- 
bution of results is regular, which indicates that our results 
are not biased by a particular object. 

V99 also investigated the deviations from spherical sym- 
metry by measuring the value of /3 in six sectors in the ra- 
dial range r > 0.3riso, and concluded that the assumption 
of spherical symmetry is relatively well satisfied in cluster 
outskirts, at variance with our results (see Fig. 7). However, 
when fitting a /3-model the fit is mostly driven by the shape 
of the profile in the innermost region, where the statistics 
are higher. Conversely, our method is model-independent, 
and it directly stacks the data at similar radii. For relaxed 
objects, our data also indicate little deviation from spheri- 
cal symmetry at r < 7500 > and a significant scatter is only 
observed beyond 7-500 > so it is probable that these devi- 
ations would not be reflected in the /3-model fit. For in- 
stance, the case of A2029 is striking. While, in agreement 
with V99, we find little azimuthal variations in /3 outer, we 
observe a high level of scatter in this object beyond r 5 oo, 
which is explained by the presence of a possible filament 
connecting A2029 to its neighbor A2033 in the north (see 
Gastaldello et al, (2010) and Appendix B). Moreover, V99 
deliberately excluded a number of systems with obviously 
disturbed morphologies, such as A3558 and A3266, which 
we included in our sample. Therefore, our results do not 
contradict those of V99. 


5. Comparison with numerical simulations 

In this section, we compare our observational re- 
sults with three different sets of numerical simulations 
(Roncarelli et al. 2006; Nagai et al. 2007b; Vazza et al. 
2010 ). We analyze the results of a composite set of cosmo- 
logical runs, obtained by the different authors with slightly 
different cosmological and numerical setups. In addition, 
the preliminary data reduction was made on each dataset 
following independent post-processing techniques, aimed at 
assessing the role of gas clumping on the comparison be- 
tween simulated mock and real X-ray observations. Our 
aim in this project is to test the most general and converg- 
ing findings of such different runs against our observations 
with ROSAT/ PSPC. 

5.1. Simulations 

5.1.1. ENZD 

We use a sample of 20 simulated clusters from the high- 
resolution and non-radiative (NR) resimulations of massive 
systems presented in Vazza et al. (2010). In this set of 
simulations, adaptive mesh refinement in the ENZ0 1.5 code 
(Norman et al. 2007) has been tailored to achieve high 
resolution in the innermost regions of clusters (following 
the increase in gas and DM over density), and also in the 
outermost cluster regions, following the sharp fluctuations 
of the velocity field, associated with shocks and turbulent 
motions in the ICM. For a detailed presentation of the 
statistical properties of the thermal gas (and of turbulent 
motions) in these simulated systems we refer the reader to 
Vazza et al. (2010, 2011a). 


5.1.2. ART 

We analyze a sample of ten simulated clusters with 
Tx > 2.5 keV from the sample presented in Nagai et al. 
(2007a, b). These simulations are performed using the adap- 
tive refinement tree (ART) N-body-hgas-dynamics code 
(Kravtsov 1999; Kravtsov et al. 2002), which is a Eulerian 
code that uses adaptive refinement to achieve high spatial 
resolution (a few kpc) in self-consistent cosmological simu- 
lations, To assess the impact of cluster physics on the ICM 
properties, we compared two sets of clusters simulated with 
the same initial conditions but with different prescription of 
gas physics. In the first set, we performed hydrodynamical 
cluster simulations without gas cooling and star formation. 
We refer to this set of clusters as NR clusters. In the second 
set, we turned on the physics of galaxy formation, such as 
metallicity-dependent radiative cooling, star formation, su- 
pernova feedback, and a uniform UV background. We refer 
to this set of clusters as cooling+star formation (CSF) clus- 
ters. For detailed descriptions of the gas physics and mock 
X-ray images we refer the reader to Nagai et al. (2007a, b). 

Following Nagai k Lau (2011), we also computed the 
dumping-corrected gas density profiles of X-ray emitting 
gas with T > 10 6 K for comparisons with X-ray observa- 
tions. Indeed, the formation of dense clumps increases the 
emissivity of the gas, which leads to an overestimation of 
the measured gas density when the assumption of constant 
density in each shell is made. For these profiles, we com- 
puted the average squared density from the simulations in 
each radial bin and took the square root of the total to 
mimic the reconstruction of density profiles from real data 
(see Nagai k Lau 2011, for details). 

5.1.3. GADGET 

This set includes four massive halos simulated with the 
GADGET-2 Tree-SPH code (Springel 2005), with M 200 > 
10 15 Mq (for a detailed description see Roncarelli et al. 
2006, and references therein). Each object was simulated 
following two different physical prescriptions: a NR run (re- 
ferred to as ovisc in Roncarelli et al. 2006) and a run in- 
cluding cooling, star formation, and supernovae feedback 
(CSF). 

To eliminate the dense clumps that dominate the den- 
sity and surface brightness in the outskirts, when comput- 
ing the profiles for every radial bin, we excise the one per 
cent of the volume that corresponds to the densest SPH 
particles. This empirical method mimics the procedure of 
masking bright isolated regions from the analysis of ob- 
served clusters. 

5.2. Comparison of gas density profiles 

We compared the simulations with our observed mean 
ROSAT density profile (see Fig. 5 and Table 1). We present 
the detailed comparison in Fig. 8, with the NR simulations 
(left hand panel) and with the CSF simulations (right). 
From the figures, we find relatively good agreement be- 
tween all the different sets of simulations, especially beyond 
~ 0.77-200- The NR GADGET run has a lower normalization 
than the corresponding grid codes, because in GADGET the 
fraction of baryons virializing into clusters is less than the 
cosmic value (~ 78% of the cosmic baryon fraction), while 
grid codes predict a baryon fraction in clusters very close 
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Fig. 8. Comparison between the mean ROSAT density profile for our sample and the different sets of numerical simula- 
tions. The shaded area indicates the data and la scatter as shown in Fig. 5. The bottom panels show the ratio between 
simulations and data as a function of radius. Left: Comparison with NR simulations. The dotted red curve represents the 
ENZQ profile (Vazza et al. 2010), the solid green curve shows the ART simulations (Nagai et al. 2007b), and the dashed 
blue curve is the GADGET profile (Roncarelli et al. 2006). Right: Same with CSF simulations. The dashed blue line shows 
the GADGET simulations, while the green curves show the ART profiles, for the total density (solid) and corrected for 
clumping (dotted, Nagai k Lau 2011). 


to the cosmic value. In general, we see that the predicted 
density profiles are too steep compared to the data. We 
note that NR runs predict steeper profiles than the runs in- 
cluding cooling, star formation, and feedback effects. CSF 
profiles also have lower normalizations, since radiative cool- 
ing transforms a fraction of the gas into stars. The profile 
including the effects of clumping shows the best agreement 
with the data. 

To quantify this effect, we fitted the various profiles in 
three different radial ranges (0.2 — 0.4r2och 0*4 — 0.65r2oo, 
and 0.65 — 1.2r*2oo)* In the inner regions, the effects of addi- 
tional physics are expected to be important, thus highlight- 
ing the differences between NR and CSF runs. The radial 
range 0.4 — 0.65r2oo(^ 0.6 — lrsoo) is a good range for com- 
paring with the data, since the effects of radiative cooling 
should be small, and data from several different satellites 
are available for cross-check. On the observational side, the 
density profiles in this radial range are well-fitted by the ft- 
model (see Eq. 4), and several independent works converge 
to the canonical value of ft ~ 0.7 (e.g., Mohr et al. 1999; 
Ettori k Fabian 1999; Vikhlinin et al. 1999; Croston et al. 
2008; Ettori & Balestra 2009; Eckert et al. 2011b). As a 
benchmark, we computed the values of ft for our average 
density profile and the various sets of simulations, fixing the 
core radius to 0.12r2oo (e.g., Mohr et al. 1999). The results 
of this analysis are shown in Table 2. The fits to the ob- 
servational data were performed on the emission-measure 
profiles (see Sect. 4.2) to take advantage of the larger num- 
ber of bins and minimize the uncertainties linked to the 
deprojection procedure. 

These numbers confirm the visual impression that the 
simulated gas density profiles are steeper than the observed 
ones. In the 0.4 — 0.65r , 2oo range, while all our datasets 
converge to a ft value very close to the canonical value, 
all the simulations lead to significantly steeper gas profiles, 


Table 2. Values of the ft parameter 

(Cavaliere k Fusco-Femiano 1976) in several radial 
ranges for the average ROSAT profiles and the various 
sets of simulations. 


Data set 

00.2-0.4 

00.4-0.65 

00.65 — 1.2 

Data, total 

0.661 ± 0.002 

0.710 ± 0.009 

0.890 ±0.026 

Data, CC 

0.700 ± 0.004 

0.699 ± 0.016 

1.002 ±0.057 

Data, NCC 

0.635 ± 0.003 

0.723 ± 0.011 

0.853 ±0.029 

ENZQ 

0.744 

0.945 

0.952 

ART, NR 

0.801 

0.956 

0.983 

ART, CSF 

0.808 

0.842 

1.005 

ART, NR, cl 
ART, CSF, cl 

0.701 

0.824 

0.854 

0.803 

0.718 

0.902 

GADGET, NR 

0.856 

0.857 

0.971 

GADGET, CSF 

0.756 

0.864 

0.944 


Note: The core radius was fixed to 0.12r2oo in all cases. The 
subscript cl indicates the profiles corrected for the effect of 
clumping using the method described in Nagai k Lau (2011). 


with ft values higher than 0.85, with the exception of the 
ART profile that includes CSF and clumping. Therefore, we 
can see that at this level of precision the effects of additional 
physics cannot be neglected, even in regions well outside of 
the cluster core. 

The results presented in Table 2 also highlight the dif- 
ferences between NR and CSF runs. Inside r$ oo, the sim- 
ulations including additional physics lead to flatter den- 
sity profiles compared to the NR runs. In this case, gas 
cooling converts a fraction of the X-ray emitting gas into 
stars. Since the cooling efficiency decreases with radius, 
more gas disappears from the X-ray range in the central 
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Fig, 9. Left: Comparison between the average observed azimuthal scatter profile from Fig. 7 (black) and the scatter in 
the simulations for the ENZO runs (red), for the total scatter (solid line) and when filtering out the 1% most-luminous 
cells (dashed curve). The cyan (NR) and magenta (CSF) curves represent the scatter in the ART simulations. Right: Same 
for the CC (red) and NCC (blue) observed profiles, compared to the 1%-filtered ENZO profiles for the morphologically 
relaxed (red) and disturbed (blue) simulated clusters. 


regions, which results in flatter density profiles and lower 
normalizations. We note, however, that this effect is prob- 
ably overestimated in the CSF simulations. Indeed, it is 
well-lmown that these simulations predict a stellar fraction 
that is well above the observed value (e.g., Kravtsov et al. 
2005; Borgani Sz Kravtsov 2009). This effect is particularly 
strong in the ART CSF simulation, for which nearly one 
third of the gas is converted into stars. Beyond 7 * 500 , there 
is little difference between NR and CSF runs; i.e., the effects 
of additional physics are not important. At large radii, the 
effect of gas clumping (Nagai &; Lau 2011) dominates and 
flattens the observed profiles. As we can see in Table 2 and 
in the right hand panel of Fig, 8, the ART profile including 
both additional physics and a post-processing treatment of 
clumping reproduces the behavior of the data more closely, 
even though it is still slightly too steep. 


5.3. Azimuthal scatter 

A study of the azimuthal scatter in the radial profiles of 
density, temperature, entropy and X-ray brightness of sim- 
ulated ENZO clusters has been presented in Vazza et al. 
(2011b). In this case, we differ from the analysis reported 
there by computing the azimuthal scatter from more an- 
gular sectors, N=12, than for N=2, 4, and 8 explored 
in Vazza et al. (2011b). In the simulations, several dense 
clumps are present, which may bias the predicted scatter. 
To overcome this problem, we computed the scatter of the 
simulated clusters both for the total gas distribution and 
by filtering out the 1% most X-ray luminous cells, as in 
Roncarelli et ai. (2006), which removes a large fraction of 
the clumps. 

We also performed a similar analysis on the set of ART 
simulations, both for the NR and CSF rims. In this case, we 
analyzed mock X-ray images using the same method as the 
observational data (see Sect. 3.4), and applied our point- 


source detection algorithm to remove the most prominent 
clumps. 3 

In Fig. 9, we show the measured scatter profile from 
Fig. 7, together with the scatter profiles of X-ray bright- 
ness from ENZO and ART simulations. Interestingly, we note 
that NR runs (red and cyan) overestimate the observed az- 
imuthal scatter, while CSF simulations underestimate it. 
In the latter case, radiative cooling lowers the entropy of 
the gas, which makes it sink into the cluster’s potential 
well. This effect produces more spherical X-ray morpholo- 
gies, thus lowering the azimuthal scatter. Conversely, in NR 
runs, the effects of dynamics are more important, which cre- 
ate more substructures and increases the azimuthal scatter. 

Interestingly, the profile that best reproduces the data 
is the ENZO profile for which the 1% most-luminous pix- 
els were filtered out. This may indicate that some clumps 
are indeed present in the observations, but were detected 
as point sources and were masked for the analysis. We note 
that, even if in this case the azimuthal scatter from NR sim- 
ulation runs is in good agreement with the ROSAT data, 
the absolute profiles of density are too steep compared to 
observations (see the left hand panel of Fig. 8). However, 
our definition of the azimuthal scatter (Eq. 2) is normal- 
ized to the absolute value of the profile at each radii, which 
makes it a rather robust proxy of cluster asymmetries on 
large ~ Mpc scales. 

In the right hand panel of Fig. 9, we also show the aver- 
age radial trends of the azimuthal scatter for the projected 
X-ray emission from the ENZO clusters after dividing the 
dataset into 11 CC-like and 9 NCC-like objects, compared 
to the observed scatter profiles for the CC and NCC classes 
from Fig. 7. This division is of course only qualitative, since 
no radiative cooling is modeled in these runs. However, our 


3 Because of hov r few objects are considered, we ignored the 
GADGET simulations for this analysis. For a comparison be- 
tween GADGET and ENZO scatter profiles, we refer the reader to 
Vazza et al. (2011b). 
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sample can be divided into classes that are quite similar 
to observed CC and NCC properties , based on the anal- 
ysis of the power ratios P3/P0 and of the centroid shift 
w, evaluated within 7*500 as in Cassano et al. (2010). We 
classify as NCC-like systems those for which the values of 
P 3 /P 0 > HT 7 and w > 0.02 were found in at least two of 
the three projected maps along the coordinate axes, or as 
CC-like otherwise, identical to what was done for the same 
sample in Vazza et al. (2011a). 

In this figure, we can clearly see that the radial trend 
of the difference between the two populations disagrees. 
While in simulations the two trends detach as we move 
farther out in the cluster atmospheres, in the observed pro- 
files the most prominent differences are found in the range 
0.2 < r/r2oo < 0.8. I11 the CC case, we find better qualita- 
tive agreement in the outskirts than in the central regions. 
This is not surprising, given that radiative cooling and en- 
ergy feedback from central AGNs are missing in these runs. 
Indeed, as we can see in the left hand panel of Fig. 9, radia- 
tive cooling has a strong impact on the general morphology 
of clusters (Fang et al. 2009; Lau et al. 2011). On the other 
hand, the simulated disturbed systems have a larger scatter 
in the outskirts than the observed NCC clusters. However, 
we observe large differences in the scatter between the var- 
ious NCC profiles, such that the result may be affected 
by small- number statistics. In any case, since the selection 
criteria are very different, we do not expect a one-to-one 
correlation between the various classes. 


6. Discussion 

6.1. Observational results 

In agreement with earlier works using ROSAT (V99, N05) 
and Chandra (Ettori & Balestra 2009), but at variance 
with some recent results from Suzaku (Bautz et al. 2009; 
Simionescu et al. 2011; George et al. 2009) and XMM- 
Newton (Urban et al. 2011), our analysis reveals that on 
average the slope of the density profiles steepens beyond 
7*5 00 (see Table 2). This result indicates that the latter re- 
sults may have been performed along preferential directions 
connected with the large-scale structure (e.g., in the di- 
rection of filaments). Indeed, the narrow FOV of Suzaku 
onlv allowed sparse coverage of the outskirts of nearby 
clusters, so that these measurements might be the result 
of azimuthal variations. In the case of A 1795, Bautz et al. 
(2009) detected a significant signal only in the northern di- 
rection, while the Perseus result (Simionescu et al. 2011) 
was obtained along two narrow arms, covering less than 
10% of the cluster’s extent at 7*200- Moreover, using sev- 
eral offset ROSAT /PSPC pointings of the Perseus cluster, 
Ettori et al. (1998) observed clear azimuthal variations in 
the density and gas fraction. Therefore, it is likely that the 
aforementioned measurements are not representative of the 
cluster as a whole. This picture is supported by our analysis 
of azimuthal variations in cluster outskirts, which suggests 
that even CC clusters exhibit significant departures from 
spherical symmetry around r 2 oo- Consequently, a full az- 
imuthal coverage is required to study the global behavior 
of cluster outer regions. 

An important result of this work is the systematic dif- 
ference between CC and NCC cluster populations observed 
beyond ~ 0.3r 2 oo (see Fig. 4). As explained in Sect. 4, this 
effect seems to be an intrinsic difference between the two 


classes, since it is does not correspond to a biased distribu- 
tion of our sample in temperature or redshift. Our scaled 
gas mass profiles provide a natural explanation for this re- 
sult (see Fig. 6). Indeed, when the appropriate scaling is 
applied, the steeper density profiles of CCs in the outskirts 
compensate exactly for the excess density in the central re- 
gions, such that clusters with the same virial mass have the 
same gas mass enclosed within 7*200, albeit distributed in 
a different way for relaxed and disturbed objects. This re- 
sult was expected in the old cooling-fiow scenario (Fabian 
1994), in which radiative cooling causes the gas to flow in- 
wards and accumulate in the central regions. While in the 
central regions AGN feedback prevents the gas from cooling 
below a certain level (e.g., McNamara & Nulsen 2007), the 
entropy injected by the central AGN is not strong enough 
to balance the flow in the outer regions of clusters, which 
explains the steep density profiles seen in Fig. 4. Conversely, 
merging events are capable of injecting a very large amount 
of energy in the ICM, which results in an efficient redistri- 
bution of the gas between the core and the outer regions 
and creates the flatter density profiles measured for NCC 
clusters. 

We also determined the typical scatter in surface- 
brightness as a function of radius (see Fig. 7) and split the 
data into the CC and NCC classes. In the central regions, 
we observe a systematic difference between CC and NCC 
clusters, with NCC clusters showing a higher level of scat- 
ter than CC. This result is easily explained by the larger 
number of substructures generally observed in NCC clusters 
(e.g., Sanderson et al. 2009). For CC clusters, we measure 
a scatter of 20% - 30% below 0.5r 2 oo> which corresponds 
to small variations (~ 10%) in gas density. This indicates 
that the azimuthal scatter in the inner regions (7* < 0.5r 2 oo) 
can be used to estimate the X-ray state of clusters, as sug- 
gested by Vazza et al. (2011b). Conversely, the scatter of 
CC profiles increases in cluster outskirts, and there is no 
observed difference between the two classes. Interestingly, 
we note that for CC clusters the turnover in Fig. 7 occurs 
around 7*500 , which coincides with the radius beyond which 
large scale infall motions and filamentary accretions are 
generallv non-negligible (e.g., Evrard et al. 1996). Inside 
rsoo, the gas is virialized in the cluster’s potential well, 
and it shows only little deviations from spherical symme- 
try. Beyond 7*500, accretion processes are important, and 
the gas is located mostly along preferential directions (i.e., 
filaments). As a result, the distribution of the gas becomes 
strongly anisotropic, even for clusters that exhibit a relaxed 
morphology in their inner regions. 

6.2. Comparison with simulations 

Comparing our density profiles with numerical simulations, 
we find that all NR simulations predict very steep profiles 
already starting from ~ 0.2r 2O o, with values of the ft param- 
eter greater than 0.85 in the 0.4 — 0.65r 2 oc range (see the 
left hand panel of Fig. 8 and Table 2). This indicates that 
including non-gravitational effects is needed to reproduce 
the observed slope, even well outside of cluster cores. The 
runs including additional physics are in better qualitative 
agreement with the observations (see the right hand panel 
of Fig. 8), although their gas fraction is too low because 
of overcooling (~ 10% compared to ~ 15%). However, it 
seems unlikely that star formation and galactic winds (as 
in the CSF runs explored here) are the only feedback mech- 
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anisms needed to reproduce observed clusters. Indeed, sim- 
ple feedback models still face severe problems in matching 
the properties of the stellar components inside galaxy clus- 
ters, as well as the properties of galaxies within them (e.g., 
Borgani & Kravtsov 2009, for a recent review). 

As illustrated in Table 2, gas clumping may also play 
a role in reconciling simulations with observations. Indeed, 
if an important fraction of the gas in cluster outskirts is 
in the form of dense gas clumps, as suggested in simula- 
tions (Nagai &; Lau 2011), the emissivity of the gas would 
be significantly increased, thus leading to an overestima- 
tion of the gas density when the assumption of constant 
density in each shell is made. Our results show that the 
treatment of gas clumping slightly improves the agreement 
between data and simulations (see the right hand panel of 
Fig. 8). In addition, gas clumping also provides an alterna- 
tive interpretation for our observed difference between the 
CC and NCC populations beyond 0.3r2oo- Indeed, simula- 
tions predict a larger clumping factor in unrelaxed clusters 
compared to relaxed systems for the same average density, 
which would result in a higher observed densitv in the for- 
mer. At the moment, it is not clear whether this difference 
is caused by gas redistribution or clumping, or if both of 
these effects play a role to some extent. 

On the other hand, we find that numerical simulations 
can reproduce qualitatively the observed azimuthal scatter 
in the galaxy cluster gas density profiles (see Fig. 9), al- 
though they fail to reproduce the trends observed for the 
CC and NCC populations separately. Interestingly, we find 
that the observed azimuthal scatter is reproduced with rea- 
sonable accuracy when the 1% most luminous clumps are 
filtered out, whereas the NR simulations with no filtering 
overestimate the observed level of azimuthal scatter at all 
radii. Two possible interpretations can be put forward to 
interpret this result. Observationally, it is possible that the 
dense clumps were detected as point sources and were fil- 
tered out of our observations. If this is the case, long ex- 
posures with high- resolution X-ray telescopes ( Chandra or 
XMM-Newton) should allow us to characterize the point 
sources and distinguish between dense clumps and back- 
ground AGN, possibly unveiling the population of accreting 
clumps in cluster outskirts. Conversely, if such observations 
do not confirm the existence of the clumps, it would imply 
that NR simulations significantly overestimate the amount 
of clumping in cluster outskirts, which would weaken the 
case for the interpretation recently put forward to explain 
the flattening of the entropy profiles observed in a few cases 
(Simionescu et al. 2011; Urban et al. 2011). 

As shown in Fig. 9, radiative cooling may also help rec- 
oncile the NR simulations with the data. Indeed, radiative 
cooling lowers the entropy of the gas and makes it sink 
into the potential well, which produces clusters with more 
spherical morphologies (Lau et al. 2011) and thus reduces 
the azimuthal scatter. Since we know that this effect is over- 
estimated in our CSF simulations, radiative cooling likely 
reduces the azimuthal scatter with respect to NR simula- 
tions, although not as much as what is predicted here. This 
effect may also explain why NR simulations fail to repro- 
duce the average scatter profiles of CC clusters (see the 
right hand panel of Fig. 9). 

Alternatively, AGN feedback may be an important in- 
gredient that is rarely taken into account in numerical sim- 
ulations. Recently, Pratt et al. (2010) observed an anti- 
correlation between entropy and gas fraction, such that 


multiplying cluster entropy profiles by the local gas frac- 
tion allows recovery of the entropy profiles predicted from 
adiabatic compression; i.e., the excess entropy observed in 
cluster cores is balanced by a lower gas fraction, and the to- 
tal entropy follows the predictions of gravitational collapse. 
Mathews & Guo (2011) interpret this result in terms of the 
total feedback energv injected in the ICM through various 
giant AGN outbursts, which they estimate to be as large 
as 10 63 ergs. In this scenario, feedback mechanisms are pre- 
venting the gas from collapsing into the potential well, caus- 
ing a deficit of baryons in the inner regions of clusters, hence 
flattening the observed density profiles. Moreover, it is well 
known that this mechanism also takes place on group and 
galaxy scales, leading to shallower density profiles in the 
accreting clumps. As a result, the gas distribution in clus- 
ter outskirts would be more homogeneous than predicted in 
NR simulations, in agreement with our observed azimuthal 
scatter profiles. Therefore, although its implementation into 
numerical simulations is challenging (Sijacki et al. 2008), 
AGN feedback could be an important effect for reconciling 
simulations with observations. A more complex picture of 
the ICM, possibly including the detailed treatment of mag- 
netic fields, cosmic rays, and thermal conductions (and of 
the instabilities arising from these ingredients), would still 
represent a challenge for current cosmological simulations. 

7. Conclusion 

In this paper, we have presented our analysis of a sam- 
ple of local ( z = 0.04 - 0.2) clusters with ROSAT /PSPC, 
focusing on the properties of the gas in cluster outskirts. 
We then compared our observational results with numeri- 
cal simulations (Roncarelli et al. 2006; Nagai & Lau 2011; 
Vazza et al. 2011b). Our main results can be summarized 
as follows. 

We observed a general trend of steepening in the ra- 
dial profiles of emission-measure and gas density be- 
yond ~ r 50 o, in good agreement with earlier works 
from Vikhlinin et al. (1999), Neumann (2005), and 
Ettori & Balestra (2009). As a result, the shallow den- 
sity profiles observed in several clusters by Suzaku 
(Bautz et al. 2009; Simionescu et al. 2011) are proba- 
bly induced by observations in preferential directions 
(e.g., filaments) and do not reflect the typical behavior 
of cluster outer regions. 

- We found that NCC clusters have in average a higher 
density than CC systems beyond ~ 0.3r2oo, which can- 
not be easily explained by any selection effect. We in- 
terpreted this result by a different distribution of the 
gas in the two populations: the well-known density ex- 
cess in the core of CC clusters is balanced by a slightly 
steeper profile in the outskirts, which leads to the same 
gas mass enclosed within r 2 oo hi the two populations 
(see Fig. 6). Alternatively, this result could be caused 
by a larger clumping factor in disturbed objects, leading 
to an overestimate of the gas density of NCC clusters 
in the external regions. 

— We also observed that NCC systems have higher az- 
imuthal scatter than CCs in the central regions, which 
is easily explained by the more disturbed morphology 
of NCC clusters. Conversely, beyond ~ r 50 o both pop- 
ulations show a similar level of asymmetry (60-80%), 
which suggests that a significant fraction of the gas is 
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in the form of accreting material from the large-scale 
structures. 

— Comparing our ROSAT density profile with numerical 
simulations, we found that all NR numerical simula- 
tions fail to reproduce the observed shape of the den- 
sity profile, predicting density profiles that are signif- 
icantly too steep compared to the data (see Table 2 
and Fig. 8). This implies that non gravitational effects 
are important well outside the core region. The runs in- 
cluding additional physics (cooling, star formation, SN 
feedback) predict flatter profiles, although still too steep 
compared to the observations. Besides, it is well known 
that these simulations overpredict the stellar fraction in 
clusters (Borgani & Kravtsov 2009). A slightly better 
agreement is found when a treatment of the observa- 
tional effects of gas clumping is adopted (Nagai & Lau 
2011 ). 

NR simulations are able to predict the observed az- 
imuthal scatter profile with reasonable accuracy, but 
only when the 1% most luminous cells are filtered out 
(see Fig. 9). This result implies that either (i) the clumps 
are quite bright and were masked as point sources in our 
analysis pipeline, in which case offset X MM- Newton and 
Chandra observations will be able to characterize them 
spatially and spectrally, or (ii) the non-radiative simu- 
lations significantly overestimate the effects of clumping 
on the observable X-ray properties. Because of the ab- 
sence of cooling, it is however hard for these simulations 
to reproduce the observed trends of azimuthal scatter 
for the two populations (CC and NCC) separately. 

As an alternative explanation, we suggest that AGN 
feedback might be important even at large radii, and could 
help to reconcile observations and simulations. Indeed, re- 
cent works (Pratt et al. 2010; Mathews & Guo 2011) indi- 
cate that feedback mechanisms may be responsible for the 
well-known deficit of baryons in cluster cores, thus leading 
to flatter gas distributions out to large radii. Moreover, the 
existence of such mechanisms on group and galaxy scales 
could also dilute the accreting material at large radii, lead- 
ing to a smaller azimuthal scatter. 
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Tkble 3. Master table of the cluster sample. Column description: 1. Cluster name; 2. Effective exposure of the PSPC observation; 3. Redshift (from NED); 4. 
Hydrogen column density, Nh , along the line of sight (Kalberla et al. 2005); 5. Mean temperature in the 200-500 kpc radial range; 6. 7*200 from Arnaud et al. 
(2005) scaling relations, in physical units; 7. Same as 6, in apparent units; 8. Central density no (this work); 9. Central entropy Kq, from Cavagnolo et al. (2009); 
10. Reference for the temperature profile (l=Snowden et al. (2008); 2=Cavagnolo et al. (2009); 3=De Grandi & Molendi (2002)). 


Cluster 

Exposure [ks] 

z 

Nh [10 22 era 2 ] 

AX200-500 [keV] 

r 2 oo [kpc] 

^200 [arcmin] 

no [10 3 cm 3 ] 

Ko [keV cm 2 ] 

Reference 

A85 

10.065 

0.05506 

0.028 

6.3 ±-0.1 

1873 

29.17 

18.9 ± 0.25 

12.5 

1 

A119 

14.758 

0.0442 

0.037 

5.0 ±0.1 

1673 

32.04 

2.1 ± 0.34 

233.9 

2 

A133 

19.429 

0.0566 

0.0164 

4.0 ±0.09 

1494 

22.68 

14.0 ±0.18 

17.3 

1 

A401 

7.519 

0.07366 

0.0995 

7.9 ±0.15 

2077 

24.72 

5.3 ± 0.66 

166.9 

2 

A478 

23.019 

0.0881 

0.131 

6.56 ±0.08 

1883 

19.05 

18.8 ±0.19 

7.8 

1 

A644 

10.310 

0.0704 

0.0750 

7.7 ±0.1 

2054 

25.48 

9.4 ± 0.29 

132.4 

2 

A665 

37.066 

0.1819 

0.0431 

8.0 ± 0.2 

1987 

10.82 

5.6 ±0.18 

134.6 

1 

A1068 

10.822 

0.1375 

0.0173 

4.9 ±0.17 

1587 

10.89 

15.0 ±0.24 

9.1 

1 

A1651 

7.630 

0.084945 

0.0156 

6.7 ± 0.2 

1913 

20.00 

8.8 ± 0.50 

89.5 

2 

A1689 

14.291 

0.1832 

0.0186 

9.2 ± 0.2 

2126 

11.51 

13.8 ±0.22 

78.4 

1 

A1795 

35.494 

0.06248 

0.0121 

6.02 ± 0.08 

1828 

25.31 

20.1 ±0.12 

19.0 

1 

A1991 

21.956 

0.0586 

0.0248 

2.4 ± 0.1 

1064 

15.64 

16.1 ± 0.22 

1.5 

1 

A2029 

13.089 

0.07728 

0.0323 

7.7 ± 0.2 

2054 

23.40 

20.2 ±0.20 

10.5 

1 

A2142 

19.410 

0.0909 

0.0383 

9.0 ±0.3 

2209 

21.73 

10.3 ± 0.17 

68.1 

3 

A2163 

7.267 

0.203 

0.109 

18.8 ± 1.3 

3008 

15.01 

8.2 ± 0.92 

438.0 

2 

A2204 

5.346 

0.1526 

0.0561 

8.3 ± 0.2 

2057 

12.93 

33.3 ± 0.76 

9.7 

1 

A2218 

43.179 

0.1756 

0.0266 

6.7 ±0.3 

1825 

10.22 

4.6 ± 0.10 

288.6 

1 

A2255 

13.676 

0.0806 

0.0250 

6.1 ±0.1 

1817 

19.9 

2.3 ± 0.32 

529.1 

2 

A2256 

17.000 

0.0581 

0.0418 

6.2 ±0.1 

1865 

27.63 

3.0 ± 0.47 

349.6 

1 

A2597 

7.426 

0.0852 

0.0246 

3.64 ± 0.06 

1405 

14.65 

18.0 ±0.22 

10.6 

1 

A3112 

7.829 

0.07525 

0.0137 

4.8 ±0.1 

1613 

18.82 

18.3 ±0.26 

11.4 

1 

A3158 

3.123 

0.0597 

0.0138 

5.1 ±0.1 

1681 

24.27 

3.8 ± 0.20 

166.0 

1 

A3266 

13.967 

0.0589 

0.0158 

9.2 ±0.3 

2260 

33.05 

5.3 ± 0.49 

72.5 

3 

A3558 

28.751 

0.048 

0.0402 

5.06 ± 0.05 

1687 

29.89 

7.2 ± 0.23 

126.2 

1 

A3562 

20.518 

0.049 

0.0376 

4.8 ± 0.3 

1635 

28.41 

5.7 ±0.26 

77.4 

3 

A3667 

12.462 

0.0556 

0.0452 

5.31 ±0.05 

1721 

26.56 

4.5 ± 0.36 

160.4 

2 

A4059 

5.684 

0.0475 

0.0122 

4.07 ± 0.08 

1513 

27.08 

4.7 ± 0.33 

7.1 

1 

Hydra A 

18.541 

0.0539 

0.0468 

4.0 ±0.06 

1495 

23.75 

22.1 ±0.17 

13.3 

1 

MKW 3s 

9.781 

0.045 

0.0272 

3.52 ±0.06 

1409 

26.54 

13.5 ± 0.22 

23.9 

1 

PKS 0745-191 

9.627 

0.1028 

0,405 

8.4 ±0.3 

2121 

18.70 

31.9 ± 0.45 

12.4 

1 

Triangulum 

7.343 

0.051 

0.114 

8.9 ±0.2. 

2229 

37.31 

5.9 ±0.79 

313.0 

1 
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Appendix A: Determination of azimuthal scatter 
profiles 

The azimuthal scatter (Vazza et aL 2011b) is defined as 
the relative scatter in surface brightness between various 
sectors (see Sect. 3.4), 


a 1 ^ (SBj - {SB)) 2 

' N ti ( S - B ) 


(A.1) 


In practice, computing this quantity is difficult, since the 
statistical fluctuations of the surface brightness introduce a 
contribution to the scatter that is actually dominant in the 
outer regions. To estimate the intrinsic Level of azimuthal 
scatter, we used two different complementary methods, 
which we describe in more detail here. 


Al. Subtraction of the statistical scatter 

Since the statistical fluctuations of the data also introduce 
a certain level of scatter, it must be noted that the quantity 
computed through Eq, A.l gives the sum of the statistical 
and intrinsic scatter, 

£ 2 = sL + ELt- '(a.2) 

The statistical scatter is given by the mean of the 

individual relative errors, 

E L t = ^E^j2> ( A - 3 ) 

and must be subtracted from Eq. 2 to estimate the level of 
intrinsic scatter. The validity of Eq. A. 3 for the statistical 
scatter was verified through a set of simulations of a source 
with no intrinsic scatter. 

The uncertainties in the scatter are then estimated 
through Monte Carlo simulations. Namely, the surface- 
brightness values in the N sectors are randomized, and the 
scatter is recomputed each time. This procedure is applied 
10 3 times, and the error on the scatter is defined as the 
RMS of the distribution around the mean value. 


A.2 . Maximum likelihood estimation 

To check the validity of our approach we performed an inde- 
pendent analysis of the scatter. We model the intrinsic scat- 
ter in the form of a Gaussian. We use a maximum likelihood 
algorithm (Maccacaro et al. 1988) to fit the data, where the 
free parameters are the mean and the intrinsic scatter (i.e. 
the standard deviation of the Gaussian). The methods de- 
scribed in both Sect. A.l and this appendix were applied 
to the surface brightness distribution within the annuli of 
each cluster (see Sect 3.4 for details). Intrinsic scatter pro- 
files from different objects were rebinned onto a common 
grid in units of 7*200 and stacked. In Fig. A.l we compare 
the intrinsic scatters measured with the two methods. The 
profiles are very similar, the general trend towards increas- 
ing scatter with radius is recovered with both methods. The 
only bin where a significant difference is observed is around 
0.7^200 • This comparison therefore provides a confirmation 
of our scatter analysis using two very different methods. 


CO 



Fig. A.l. Comparison between the mean azimuthal scat- 
ter profiles computed using the direct method (black, see 
Sect. A.l) and the alternative method using a maximum 
likelihood estimator (red, see Sect. A.2). 


Appendix B: Notes on individual objects 

- A 85: 

A subcluster located ~ 10' south of the cluster center is 
currently merging with the main cluster. This substruc- 
ture was masked for the analysis. 

A401: 

The cluster is connected through a filament to its neigh- 
bor A399, located ~ 35' south-west of the center of 
A401. We extracted the surface- brightness profile in a 
sector of position angle 340-250° to avoid any contami- 
nation of A399 to our measurement of the CXB. 

- A478: 

The combination of a favorable temperature/redshift 
and a good-quality ROSAT observation allows us to 
reach the highest signal-to-noise ratio in the sample at 
7*200 for this strong CC cluster. As a result, the data from 
this cluster may contribute strongly when a weighted 
mean is performed. 

- A644: 

This NCC cluster exhibits an unusual decreasing az- 
imuthal scatter profile, showing large (close to 100%) 
scatter in its central regions, but no significant scatter 
around 7*200* 

- A2029: 

A probable filament connects A2029 to A2033, located 
~ 35' north of the center of A2029. The surface- 
brightness profile was extracted in a sector with position 
angle 140-80° to measure the CXB level. 

A2142: 

Several PSPC observations of this famous cold-front 
cluster exist. For this work, we used the longest avail- 
able observation, which was pointed 16' south of the 
center of A2142. This is the only case in the sample for 
which the observation was not pointed on the target. 

A 3558 and A 3562: 

These two clusters are located in the Shapiey super- 
cluster and connected by a filament. Consequently, they 
show an unusually high azimuthal scatter in' the out- 
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skirts. The CXB level was estimated by excluding the 
direction of the filament. 

- A 3667: 

This very disturbed cluster shows the highest emission- 
measure and density in the sample beyond ~ 0 . 2 r 2 oo> 
and hence it could bias our average profiles, in particu- 
lar when computing the difference between the CC and 
NCC classes. However, removing it from the sample did 
not lead to any significant difference, either quantitative 
or qualitative. 

- A4059: 

This is the most azimuthally-symmetric cluster in the 
sample. The azimuthal scatter for this cluster is consis- 
tent with 0 at all radii. 

Hydra A: 

A tail of emission (filament?) extends out to ~ 20' 
south-east of the cluster core. This leads to a very high 
azimuthal scatter (> 100%) around r 2 oo- 

Appendix C: Mean emission- measure profiles 

In Table C.l we give the mean self- similar scaled emission- 
measure profiles for the CC and NCC classes and the whole 
sample, as shown in Fig. 4. 


Appendix D: Computing the gas fraction from 
density profiles 

The gas fraction in the observations and in the simulated 
clusters within an overdensity A can be computed directly 
from the profiles presented in Fig. 8. Indeed, by definition, 

4 

Ma - A Pcrit^A, (D.l) 

where pcrit = f^; = 9.2 x 10 -30 g cm -3 . Then, 


p _ lv *gas, A _ 

Jgas,A ~ Ma ~ 


Aperitif 


f 


p gas {r)r 2 dr (D.2) 


Making the substitution x = we find the convenient 
formula 


fgas, A — J' Pgasip^)^ dx. 


(D.3) 


17 



D. Eckert et a!,: The gas distribution in the outer regions of galaxy clusters 


Table C.l. Data of Fig. 4: mean self-similar scaled emission-measure profiles for the whole sample and for the CC and 
NCC classes, in units of cm -6 Mpc 


Rin 

-Rout 

Total 


CC 


NCC 


0 

0.02 

(1.80 ± 0.01) • 

io - 5 

(9.48 ± 0.05) 

10“ s 

(1.13 ± 0.01) 

10 -5 

0.02 

0.04 

(1.26 ±0.01)) 

10 -5 

(4.83 ± 0.02) 

10 -5 

(8.32 ± 0.06) 

•10 -6 

0.04 

0.06 

(9.63 ± 0.04) 

10 -6 

(2.28 ± 0.01) 

10“ 5 

(6.90 ± 0.04) 

10 -8 

0.06 

0.08 

(7.39 ± 0.03) 

10 -6 

(1.23 ±0.01) 

10 _s 

(5.70 ± 0.03) 

10 -8 

0.08 

0.1 

(5.45 ± 0.02) 

10 -s 

(7.72 ± 0.04) 

10“ 8 

(4.49 ± 0.02) 

10" 8 

0.1 

0.12 

(4.12 ± 0.02) 

10 -6 

(5.27 ±0.03) 

10~ 8 

(3.52 ± 0.02) 

•10- 8 

0.12 

0.14 

(3.20 ± 1.36) 

10 -6 

(3.63 ± 0.02) 

10 “• 8 

(2.91 ± 0.02) 

•10" 8 

0.14 

0.16 

(2.47 ±0.01) 

10~ 6 

(2.60 ± 0.02) 

10" 6 

(2.37 ± 0.01) 

•10- 8 

0.16 

0.18 

(1.91 ±0.01) 

10 -8 

(1.95 ±0.01) 

10“ 6 

(1.88 ± 0.01) 

10 -8 

0.16 

0.2 

(1.51 ±0.01) 

10~ 6 

(1.48 ±0.01) 

10 -8 

(1.54 ± 0.01) 

10" 8 

0.2 

0,22 

(1.23 ±0.01) 

10“ 6 

(1.19 ±0.01) 

10~ 8 

(1.26 ± 0.01) 

10 -8 

0.22 

0.24 

(1.02 ±0.01) 

10 -8 

(9.47 ± 0.09) 

1(T 7 

(1.07 ±0.01) 

10“ 8 

0.24 

0.26 

(8.40 ±0.05) 

1(T 7 

(7.61 ± 0.08) 

1(T 7 

(8.95 ± 0.07) 

■10" 7 

0.26 

0,29 

(6.91 ±0.05) 

10“ 7 

(6.09 ± 0.07) 

10~ 7 

(7.59 ± 0.06) 

■ icr 7 

0.29 

0.31 

(5.32 ±0.04) 

10“ 7 

(4.73 ± 0.06) 

10“ 7 

(5.77 ± 0.05) 

•10- 7 

0.31 

0.34 

(4.30 ±0.04) 

10“ 7 

(3.74 ± 0.06) 

10" 7 

(4.70 ± 0.05) 

• 10 -7 

0.34 

0.38 

(3.20 ±0.03) 

10 -7 

(2.77 ±0.04) 

10“ 7 

(3.60 ± 0.04) 

•10- 7 

0.36 

0.41 

(2.49 ±0.02) 

10 -7 

(2.10 ±0.04) 

10“ 7 

(2.76 ± 0.03) 

• 10 -7 

0.41 

0.45 

(1.86 ± 0.02) 

10 -7 

(1.57 ±0.03) 

10“ 7 

(2.11 ± 0.03) 

•10 -7 

0.45 

0.50 

(1.48 ± 0.02) 

10" 7 

(1.27 ±0.03) 

10“ 7 

(1.63 ±0.02) 

•10“ 7 

0.50 

0.55 

(1.07 ± 0.02) 

10“ 7 

(9.05 ± 0.24) 

10“ 8 

(1.18 ± 0.02) 

•10- 7 

0.55 

0.60 

(7.99 ±0.14) 

10 -8 

(6.82 ± 0.22) 

10~ 8 

(8.87 ±0.19) 

•10 -8 

0.6C 

0.66 

(5.73 ±0.12) 

10 -8 

(4.97 ±0.18) 

10" 8 

(6.30 ±0.16) 

•10 -8 

0.66 

0.72 

(4.28 ±0.11) 

10 -8 

(3.78 ±0.17) 

10" 8 

(4.62 ± 0.14) 

•10“ 8 

0.72 

0.79 

(3.06 ±0.11) 

10“ 8 

(2.75 ±0.18) 

10" 8 

(3.21 ±0.13) 

•10 -8 

0.79 

0.87 

(2.23 ±0.10) 

10 -8 

(1.77 ±0,16) 

10“ 8 

(2.51 ±0.13) 

•10" 8 

0.87 

0.95 

(1.35 ± 0.09) 

10 -8 

(8.57 ± 1.49) 

10 -8 

(1.63 ±0.11) 

• 10 -8 

0.95 

1.05 

(7.77 ± 0.85) 

10“ 8 

(5.85 ± 1.40) 

10" 9 

(8.88 ± 1.07) 

•10 -9 

1.05 

1.15 

(5.32 ± 0.80) 

10 -8 

(4.19 ± 1.35) 

10“ 8 

(5.92 ± 0.99) 

•10 -9 

1.15 

1,26 

(4.74 ± 0.81) 

10 -8 

(3.75 ± 1.40) 

10“ 9 

(5.24 ± 0.97) 

• 10~ 9 
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